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РЕЗЮМЕ 


Исследованы особенности распределения численности цестод Ditestolepis diaphana (Colodkowsky 1906) и 
нематод Heligmosomum mixtum (Schulz 1952), паразитирующих в популяциях обыкновенной бурозубки Sorex 
araneus L. 1758 и рыжей полевки Myodes glareolus (Schreber 1780) соответственно. Показано, что числен- 
ность паразитов в большинстве случаев моделируется негативно-биномиальным распределением (НБР). 
Выявленные закономерности позволяют судить об устойчивом характере взаимоотношений в системе 
«паразит-хозяин». Предпринята попытка интерпретации параметров HBP, характеризующих особенности 
варьирования средней численности и дисперсии численности паразитов, связанные, по-видимому, с разли- 


чиями в восприимчивости хозяев. 
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ABSTRACT 


Abundance distribution patterns of the cestode Ditestolepis diaphana (Colodkowsky, 1906) in the common shrew 
Sorex araneus L., 1758 and the nematode Heligmosomum mixtum (Schulz, 1952) in the bank vole Myodes glareolus 
(Schreber 1780) populations were investigated. We prove that in most cases parasite abundances are modeled by 
negative binomial distribution (NBD). The patterns revealed for the parasite species studied evidence relations 
in the host-parasite system to be quite steady. An attempt was made to interpret NBD parameters describing 
variations of mean abundance and variance of the parasite abundances, which are apparently due to differences in 


host susceptibility. 
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ВВЕДЕНИЕ 


В настоящее время традиционно считается, 
что заражение паразитами имеет агрегированный 
характер, а негативно-биномиальное распреде- 
ление (HBP) – его наиболее адекватная модель 
(Crofton 1971a, b; Бреев 1972; Anderson and May 
1978; Haukisalmi 1986; Behnke et al. 1999; Ribas and 
Casanova 2005; Stear et al. 2007; Brunner and Ostfeld 
2008). Механизмы, формирующие перерассеян- 
ное (агрегированное) распределение паразитов, 
весьма разнообразны, но ведущим можно назвать 
различия в восприимчивости хозяев к заражению, 
или вариабельность дозы заражения (Anderson et 
а|. 1978). Указанный тип распределения паразита 
характеризуется тем, что большая часть хозяев в 
популяции инвазирована слабо, лишь единичные 


особи заражены значительно. Если принять за Х 
случайную величину, равную числу паразитов на 
конкретном хозяине, то отрицательное биноми- 
альное распределение имеет следующий вид: 


P{x 2i] = CL, ,p'(01- py, 

где Р обозначает вероятность того, что X = ij; 
C — число сочетаний из ik-1 по i; i = 0, 1, 2, ...; k, 
p – параметры HBP (0 « p « 1, k 70). 

Ha основании анализа встречаемости раз- 
личных видов паразитов нами (Павлов и Иешко 
1986) была построена математическая модель 
паразито-хозяинных отношений, опирающаяся на 
НБР которое определяется двумя параметрами — 
ри k. Параметр А наиболее широко обсуждается 
в современной паразитологической науке как об- 
ратный показатель агрегированности: иными сло- 
вами, чем меньше Ё, тем выше агрегированность, 
и наоборот. При малых значениях k в популяции 
хозяев доминируют особи с высокой устойчи- 
востью к заражению, в силу чего большая часть 
хозяев слабо заражена, и лишь единичные особи 
имеют высокие значения интенсивности инвазии. 
Другой же параметр НБР (р) практически не об- 
суждается. Его значения, характеризующие осо- 
бенности выживания паразитов, определяются 
отношением средней к дисперсии численности 
паразитов (p = М/О). 

В процессе эволюции паразитарных отноше- 
ний возникает некое равновесное состояние между 
популяциями паразита и хозяина, позволяющее 
паразитарной системе существовать неопределен- 
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но долгое время (Контримавичус 1982). Согласно 
Брадлей (Bradley 1972) одним из основных ме- 
ханизмов такого равновесия является регуляция 
популяций паразитов, которая осуществляется 
1) прямым влиянием факторов внешней среды, 
2) смертностью хозяев (естественной и обуслов- 
ленной воздействием паразитов) и 3) посредством 
иммунных реакций хозяина, ограничивающих 
выживаемость и плодовитость паразитов. 

Исходя из этих представлений, анализ из- 
менчивости значений интенсивности заражения 
паразитами в исследуемой популяции хозяев 
должен строиться не только на анализе средних 
значений (М), но и на оценке варьирования по- 
казателей дисперсии численности паразитов (D), 
что в какой-то мере отражает результат склады- 
вающихся регуляторных отношений в системе 
«паразит-хозяин». В этом случае параметры 
распределения численности паразитов (при НБР 
это — ри k) могут иметь биологическую интер- 
претацию, характеризуя тем самым различные 
аспекты паразито-хозяинных отношений. 

В настоящей работе делается попытка оценить 
некоторые особенности варьирования параметров 
НБР паразитов мелких млекопитающих. 


МАТЕРИАЛ И МЕТОДИКА 


Материалом для исследований послужили 
сборы паразитов мелких млекопитающих, выпол- 
ненные в 1994—2006 гг. в окрестностях Гомсель- 
гского научного стационара Института биологии 
КарНЦ РАН (62°04'N, 33*55'E), расположенного 
в среднетаежной подзоне Карелии. Основные 
исследования проводили в июне, июле и августе; 
дополнительно в отдельные годы материал соби- 
рали в мае и октябре. 

В качестве объекта изучения выбраны две па- 
разитарные системы: цестода Ditestolepis diaphana 
(Colodkowsky, 1906) — обыкновенная бурозубка 
Sorex атапеиз L., 1758 и нематода Heligmosomum 
mixtum (Schulz ,1952) — рыжая полевка Myodes 
glareolus Schreber, 1780. Для анализа основных 
параметров распределения численности были 
использованы результаты паразитологических 
вскрытий 254 экз. обыкновенной бурозубки и 
382 экз. рыжей полевки, отловленных в разные 
годы на трех мониторинговых площадках. Сбор 
материала проводили методом полного парази- 
тологического вскрытия (Аниканова и др. 2007). 


Особенности динамики паразитов 


Учитывали общее число паразитов, найденных в 
кишечнике исследованных зверьков. 

Для количественной характеристики зара- 
женности мелких млекопитающих паразитами 
использовали следующие показатели: P, % — экс- 
тенсивность инвазии (встречаемость или процент 
заражения); М, экз. — индекс обилия (средняя 
численность паразитов); дисперсия средней 
численности (D); оценки параметров HDP k и p. 
Обработку количественных данных (расчеты Р, 
М, D и *) и анализ типа распределения прово- 
дили с использованием программы Quantitative 
Parasitology (ОР) (Rozsa et al. 2000). C ee помо- 
щью проверялась гипотеза о согласовании эмпи- 
рических данных с НБР Достоверность различий 
наблюдаемых и теоретических частот распределе- 
ния принималась при Р < 0.05. Параметр р НБР 
рассчитывался по формуле: p = M/D. 
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РЕЗУЛЬГАТЫ И ОБСУЖДЕНИЕ 


D. diaphana — специфичный паразит землеро- 
ек. Жизненный цикл паразита, кроме дефинитив- 
ного хозяина, включает одного промежуточного, 
роль которого играют жуки Geotrupes stercorosus, 
Tachinus раШрез, Silpha obscura, Panorpa communis 
(larvae) u Litobius sp. (Kisielewska 1961; Vaucher 
1971; Jourdane 1975). 

Изучение индивидуальных особенностей aapa- 
жения бурозубок цестодой D. diaphana показало, 
что независимо от типа биотопа и в разные годы 
распределение численности паразита моделиру- 
ется НБР. Показатели численности и параметры 
НБР k и p представлены в Табл. 1. 

В популяции хозяина зараженность отдельных 
особей может варьировать в широких пределах. 
При агрегированном распределении значитель- 


Таблица 1. Показатели заражения и параметры НБ распределения численности цестоды Ditestolepis diaphana в популяции обык- 


новенной бурозубки. 


Table 1. Infection indices and parameters of the negative binomial distribution (NBD) of cestode Ditestolepis diaphana numbers in the 


common shrew population. 








США Li pisei en р M D k p а 
Biotope/Year Ind. investigated у with NBD) 

2 / 1997 25 36 10.96 1574.6 0.08 0.01 

3 / 1997 29 31 3.00 72.4 0.11 0.04 

1/ 1998 16 25 1.19 12.3 0.12 0.10 

З / 1998 21 43 6.52 156.9 0.15 0.04 

1/ 1999 12 75 25.75 828.0 0.37 0.03 

2 / 1999 14 50 744 263.0 0.19 0.03 

З / 1999 14 27 1.07 4.5 0.18 0.24 

2 / 2001 25 52 3.08 28.6 0.30 0.11 

3 / 2004 13 38 1.77 9.2 0.25 0.19 

2 / 2005 23 30 4.52 155.5 0.10 0.03 

3/2005 23 22 1.74 16.7 0.08 0.10 

2 / 2006 19 63 13.69 590.3 0.26 0.02 


Обозначения: P — встречаемость, %; М — средняя численность паразитов, экз.; D — дисперсия средней численности; k и p — оценки 


параметров НБР. 


Designations: P — prevalence, %; М — mean, spm.; D — variance; k and p — estimations of NBD parameters. 
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Рис. 1. Зависимость оценок параметров HBP (k, p) от встречаемости (Р, %) цестод Ditestolepis diaphana. Логарифмическая шкала. 
р: р р р Dp Ір. р 


Fig. 1. Dependence of the estimations of NBD parameters (К, p) on the prevalence (Р, %) of cestode Ditestolepis diaphana. Logarithmic 


scale. 


ная часть хозяев имеет слабое заражение, и лишь 
единичные особи заражены большим числом па- 
разитов. В силу этого зараженность бурозубок в 
популяции будет определяться параметром k HBP, 
что характеризует соотношение в популяции хо- 
зяев особей восприимчивых и не восприимчивых 
к заражению. 

Оценивая изменчивость значений параметра 
k цестод на обследованных биотопах в разные 
годы, можно отметить его значимую связь с про- 
центом заражения хозяев (Рис. 1). По мере роста 
экстенсивности инвазии (Р) агрегированность 
заражения зверьков закономерно уменьшалась. 
На рисунках приведены уравнения степенной 
регрессии и значения коэффициента детермина- 
ции (R?), достоверность которого определялась 
при Р < 0.05. Значимые значения коэффициентов 
детерминации выделены на рисунках жирным 
шрифтом. 

Для D. diaphana параметр р НБР имеет обрат- 
ную тенденцию, и с нарастанием встречаемости 
(P) его значения падают (см. Рис. 1). Если пред- 


положить, что его значения определяют прижи- 
ваемость червей, то из этого следует, что при росте 
зараженности бурозубок цестодами их приживае- 
мость падает. Значительный разброс значений р 
и невысокое значение В? свидетельствуют о об 
отсутствии достоверной связи этого параметра с 
зараженностью (Р). 

Анализ зависимости параметров р и k HBP от 
средней численности цестод (М) указывает, что 
агрегированность при росте средней интенсив- 
ности заражения меняется с меньшей скоростью, 
чем происходит снижение параметра р (Рис. 2). 
Эти особенности, на наш взгляд, связаны с тем, 
что А как показатель, который характеризует pe- 
акцию хозяина на заражение, при наблюдаемых 
значениях индекса обилия играет меньшую роль 
в регуляции численности паразитов. С другой 
стороны, обратная связь показателя р со средней 
численностью цестод и высокие значения К? 
уравнения регрессии указывает на достоверное 
и значимое снижение приживаемости цестод при 
росте интенсивности инвазии М (Рис. 2). 
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Рис. 2. Зависимость оценок параметров HBP (k, р) or средней численности (M) цестод Ditestolepis diaphana. Логарифмическая 


шкала. 


Fig. 2. Dependence of the estimations of NBD parameters (k, p) on mean abundance (M) of cestode Ditestolepis diaphana. Logarithmic scale. 


Учитывая тот факт, что распределение цестод 
в популяции бурозубки независимо от биотопа 
и года проведения исследований моделируется 
НБР можно говорить о наличии закономерной 
связи между средней численностью и ее диспер- 
сией. Согласно Тэлору (Taylor 1961) дисперсия 
численности D = а M^, гдеаи b — коэффициенты, 
определяющие взаимосвязь дисперсии и средней 
численности в сложившихся условиях взаимо- 
действия популяций паразита и хозяев. Автор 
предположил, что для устойчивых популяций жи- 
вотных значения экспоненты b близки к 2, а для 
неустойчивых — к 1. Изучение многолетней дина- 
мики численности цестоды D. diaphana показало, 
что данный параметр имеет значение, близкое к 2 
(b = 1.69), которое свидетельствует об устойчивом 
характере взаимодействия популяций цестод и 
бурозубок (Рис. 3). 

Жизненный цикл Н. mixtum до недавнего 
времени специально не изучался, но по анало- 
гии с представителями того же семейства можно 
предположить, что это — геогельминт, развитие 


которого протекает без участия промежуточных 
хозяев. Доразвитие яиц, вылупление личинок и 
линька происходят во внешней среде, во влажной 
почве или лесной подстилке, в хозяина нематоды 
попадают перорально с пищей или водой. Вид яв- 
ляется широко распространенным в Палеарктике 
паразитом лесных полевок. 

Исследования нематод, паразитирующих B ры- 
жей полевке, показали, что встречаемость парази- 
тов невысока, и только один раз была зафиксиро- 
вана 100%-ная зараженность. При обнаруженной 
невысокой в целом интенсивности инвазии и 
встречаемости нематод численность не всегда мо- 
делировалась HBP. Только для тех случаев, когда 
было обнаружено согласование с НБР, нами при- 
ведены значения параметров р и k (Табл. 2). 

Зависимость агрегированности k от встречае- 
мости (P) нематод H. mixtum в популяции полев- 
ки описывается сходным уравнением регрессии, 
что и для цестод, паразитирующих в бурозубках 
(Рис. 4). Значения показателей степени при P, или 
угол наклона регрессионных кривых, не имели до- 
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Рис. 3. Зависимость дисперсии (D) or средней численности (M) цестод Ditestolepis diaphana. Логарифмическая шкала. 


Fig. 3. Dependence of variance (D) on mean abundance (M) of cestode Ditestolepis diaphana. Logarithmic scale. 
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Рис. 4. Зависимость оценок параметров HBP (k, p) от встречаемости (P, %) нематод Heligmosomum mixtum. Логарифмическая шкала. 


Fig. 4. Dependence of the estimations of NBD parameters (k, p) on the prevalence (P, %) of nematode Heligmosomum mixtum. Logarithmic scale. 
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Таблица 2. Показатели заражения и параметры НБ распределения численности нематоды Heligmosomum mixtum в популяции 


рыжей полевки. 


Table 2. Infection indices and parameters of the NBD of nematode Heligmosomum mixtum numbers in the bank vole population. 

















Биотоп/Тод Исследовано, экз. R% M D k p (m t 
Biotope/Year) Ind. investigated with NBD) 
1/94 16 31 0.87 5.0 0.23 0.17 * 
3/95 7 100 3.43 22.9 1.19 0.15 + 
2/96 19 42 0.79 1.2 нет (по) 
2/97 19 42 0.58 0.7 нет (по) 
1/98 17 53 1.41 4.0 0.70 0.35 + 
2/98 40 40 2.47 82.4 0.17 0.03 + 
1/99 28 46 0.86 1.3 нет (по) 
2/99 12 42 0.83 2.9 0.50 0.29 + 
3/99 54 20 0.32 0.5 нет (по) 
2/00 12 42 0.83 2.9 0.50 0.29 + 
2/01 20 25 0.50 1.0 нет (no) 
2/02 23 30 0.70 1.9 0.32 0.37 + 
2/03 47 45 0.85 1.6 0.86 0.53 + 
3/03 35 31 0.60 1.1 0.49 0.53 + 
2/04 10 50 3.2 50.1 0.26 0.06 + 
2/05 12 33 1.0 3.8 0.26 0.26 + 
2/06 11 36 0.64 1.1 нет (по) 





Обозначения: P — встречаемость, %; М — средняя численность паразитов, экз.; D — дисперсия средней численности; k и р — оценки 


параметров НБР. 


Designations: P — prevalence, %; М — mean, ind.; D — variance; & and p — estimations of NBD parameters. 


стоверных различий (см. Рис. 1 и 4), что позволя- 
ет предположить сходный характер устойчивости 
хозяев к увеличению встречаемости обоих видов 
паразитов. 

Однако в отличие от цестод для H. mixtum 
вариабельность значений параметра р, как это 
видно на рис. 4, не имела выраженной связи со 
встречаемостью (P), что подтверждается недосто- 
верным значением коэффициента детерминации. 
Иными словами, при росте зараженности полевки 
нематодами их приживаемость имеет лишь незна- 
чительный тренд в сторону снижения. 


При анализе зависимости параметров k HBP 
от средней численности нематод (М) можно от- 
метить, что агрегированность № при росте интен- 
сивности заражения закономерно не меняется, 
полученные значения В? были недостоверны 
(Рис. 5). Эти особенности изменения значений k, 
по-видимому, связаны с тем, что интенсивность 
заражения неконтролируется сопротивляемостью 
хозяев, и при наблюдаемых значениях варьирова- 
ния индекса обилия хозяин не играет роли в регу- 
ляции численности как цестод, так и нематод. В то 
же время снижение параметра р является досто- 
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Рис. 5. Зависимость оценок параметров HBP (k, p) от средней численности (М) нематод Heligmosomum mixtum. Логарифмическая шкала. 


Fig. 5. Dependence of the estimations of NBD parameters (№, p) оп mean abundance (M) of nematode Heligmosomum mixtum. Logarithmic scale. 
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Рис. 6. Зависимость дисперсии (D) or средней численности (M) нематод Heligmosomum mixtum. Логарифмическая шкала. 


Fig. 6. Dependence of variance (D) оп mean abundance (M) of nematode Heligmosomum mixtum. Logarithmic scale. 


Особенности динамики паразитов 


верным, коэффициент детерминации К? равнялся 
0.609 (см. рис. 5), а обратная связь показателя р 
со средней численностью указывает на снижение 
приживаемости нематод при увеличении средней 
интенсивности заражения. Кроме того, следует 
отметить и различия коэффициентов регрессии 
(или углов наклона прямых) параметра р для 
этих паразитов (см. Рис. 2 и 5), указывающих на 
разный характер реакции паразитов при росте 
интенсивности заражения. 

Учитывая сказанное, в тех случаях, когда рас- 
пределение нематод в популяции рыжей полевки 
моделировалось НБР, мы можем характеризовать 
систему «паразит-хозяин» как устойчивую. Это 
подтверждается наличием достоверной связи 
между дисперсией и средней численностью не- 
матод (Рис. 6), где показатель степени b perpec- 
сионного уравнения равняется 2.11, т. е. имеет 
значение, близкое к 2, характеризуя тем самым 
устойчивость популяций нематод и рыжей по- 
левки. В те годы, когда не было установлено со- 
гласования распределения численности нематод в 
популяции хозяина с законом НБР (см. Табл. 2), 
коэффициент регрессии (b) дисперсии и средней 
численности имел значения, равные 0.86 (при 
В? = 0.796), что и свидетельствует о неустойчиво- 
сти паразито-хозяинных отношений популяций 
нематод в эти годы. 


ЗАКЛЮЧЕНИЕ 


В современной литературе существует много 
примеров моделирования динамики паразитарно- 
го заражения, которые основываются на использо- 
вании понятий агрегированного распределения и 
свидетельствуют о TOM, что именно гетерогенность 
хозяев по восприимчивости к заражению являет- 
ся условием формирования перерассеянного рас- 
пределения паразитов в популяции (Anderson and 
Мау 1978; May and Anderson 1978; Kennedy 1984; 
Иешко 1988; Roberts 1995). 

Следует отметить, что большинство моделей 
распределения численности паразитов было 
получено, базируясь на экспериментальных за- 
ражениях, когда удавалось контролировать усло- 
вия и дозу заражения или гетерогенность хозяев 
(возраст, пол и др.). Сравнительно небольшое 
количество работ посвящено моделированию 
заражения паразитами природных популяций 
животных. В этой связи наша работа по изучению 
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многолетней динамики численности и характера 
распределения D. diaphana в популяции обыкно- 
венной бурозубки и H. mixtum — паразитов рыжей 
полевки — показывает возможность использова- 
ния имеющихся данных для моделирования 
заражения и интерпретации параметров НБР 
Основной вывод состоит в том, что именно НБР 
является моделью, описывающей взаимодействие 
в системе «паразит-хозяин», в случае устойчиво- 
го характера взаимодействия видов. 

Существующие природные паразитарные 
системы устойчивы, и важным условием этого 
служит тот факт, что заражение паразитом не 
направлено на выработку генетически закреплен- 
ной резистентности хозяев к паразиту. Заражение 
определенным количеством паразитов является 
результатом действия двух разнонаправленных 
факторов, формально определяемых как сопро- 
тивляемость хозяина и приживаемость паразита. 
В этой связи, как показали наши исследования, 
оценки параметров k и p HBP могут быть исполь- 
зованы как показатели, отражающие динамику 
взаимодействия в системе «паразит-хозяин». 
При этом k характеризует изменения заражения 
хозяина, связанные с различиями в индивиду- 
альной устойчивости, силе иммунных реакций 
хозяев, а р определяет успех размножения или 
приживаемости паразитов. По нашему мнению, 
проведение подобных исследований имеет хоро- 
шие перспективы в развитии понимания природы 
паразитизма, его места и роли в существовании и 
динамике сообществ. 
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